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ABSTRACT 

A refinement in the satellite geopotential solution for a Goddard Earth Model 
(GEM 3) has recently been obtained. A previous satellite model, GEM 1, was 
based primarily upon satellite optical data on 25 satellites all with inclinations 
greater than 28° . The new solution includes the addition of two low inclination 
satellites, SAS at 3° and PEOLE at 15°, and is based upon 27 close earth satel- 
lites containing some 400,000 observations of electronic, laser, and optical data. 
In addition a new combination satellite/gravimetry solution (GEM 4) is derived. 
The new model includes 61 center of mass tracking station locations with data 
from GRARR, Laser, MOTS, Baker-Nunn, and NWL Tranet Doppler tracking 
sites. 

Improvement has been obtained for the zonal coefficients of the new models 
and is shown by tests on the long period perturbations of the orbits . Individual 
zonal coefficients agree very closely among different models that contain low 
inclination satellites. Tests of models with surface gravity data show that the 
GEM 3 satellite model has significantly better agreement with the gravimetry 
data than the GEM 1 satellite model, and that it also has better agreement with 
the gravitmetry data than the 1969 SAO Standard Earth II model. 
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GRAVITATIONAL FIELD MODELS GEM 3 AND 4 


I. INTRODUCTION 

Two Goddard Earth Models, GEM 1 and GEM 2, have previously been de- 
veloped and presented in a report, (1) entitled "Gravitational Field Models for 
the Earth. " The GEM 1 model was derived from satellite optical data and the 
GEM 2 model was a combined satellite and surface gravity data solution. Two 
additional solutions GEM 3 and GEM 4, employing a broader data base, are pre- 
sented here with emphasis placed on areas of improvement. The previous re- 
port should be used as a reference since a number of the techniques developed 
there are the same and are not repeated in this paper. Also reference is given 
to a number of results from the previous solutions. Table 1 presents a brief 
summary of the four Goddard Earth Models for purposes of comparison. 

1. 1 Description of GEM 3 and GEM 4 Solutions and Data Base 

A geopotential solution based upon 27 close earth satellites, including two 
low inclination satellites and some 400,000 observations of electronic, laser, 
and optical data, has recently been derived. This solution is designated as 
Goddard Earth Model, GEM 3. The spherical harmonic coefficients are com- 
plete to degree and order 12 as in the GEM 1 satellite solution which contained 
25 satellites all with inclinations greater than 28 degrees . The addition of two 
satellites PEOLE (7010901) with a 15° inclination and SAS (7010701) with a 3° 
inclination provided improved coverage and better zonal coefficients, complete 
to degree 22. The new satellite solution was combined with the surface gravity 
data for a satellite/gravimetry geopotential solution. This solution, designated 
GEM 4, has spherical harmonic coefficients complete to 16 x 16 as in the GEM 2 
solution. The geopotential coefficients of the GEM 3 and 4 solutions are listed 
in Table 2. All four models include higher degree zonal and selected satellite 
resonant coefficients extending to degree 22. 

Geocentric station coordinates were obtained for 61 tracking stations, con- 
sisting of 19 Baker-Nunn, 23 MOTS, 4 GRARR, 3 Goddard laser, and 12 NWL 
Doppler tracking sites. 

The data employed in the new GEM 3 satellite solution was based upon: 

300 weekly orbital arcs of optical data (primarily) for the 25 satellites in 

Table 3A and included 120,000 observations. 

48 weekly arcs of electronic laser, and optical data for 8 satellites, namely 

BE-B, BE-C, DI-C, DI-D, GEOS-I, GEOS-II, SAS, and PEOLE. See 

Table 3B for the distribution of data consisting of 292,000 observations. 
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100 one/two day arcs on GEOS I and II flashing light satellites, used princi- 
pally to support the MOTS station coordinates , 

The distribution of data on the 48 weekly arcs of electronic, laser, and addi- 
tional optical data is presented in Table 3B. 

The surface gravity data employed in the combination solution GEM 4 was 
the same as described in the GEM 2 solution. W Techniques and starting values 
for the solution, including the reference ellipsoid and gravity parameters, were 
the same as for GEM 1 and 2 solutions as given in the previous report. W The 
surface gravity data served as a basis for testing the satellite models and for 
comparison with other models . The data is described and tabulated in an appen- 
dix to this report. 
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II. RESULTS 


1. Zonals 


Zonal coefficients are compared among different solutions to examine the 
effect of the addition of low inclination satellite data and surface gravity data. 
Zonals are presented in Table 4 for the GEM 1, 2, 3 and 4 solutions, the SAO 
S.E. II, and the French 1971^ zonal solution. The French solution combined 
the long term zonal equations of Kozai 1969 (Kozai’s zonals being part of SAO 
S.E. II model) with corresponding equations for three low inclination satellites 
(SAS, PEOLE, DIAL). Zonal coefficient differences of each solution w.r.t. the 
French solution are presented in Table 4. The rms of the normalized coefficient 
differences scaled by 10 9 are given below: 



Solutions Including 
Effects of Low In- 
clination Satellites 

Solutions Not Including Satel- 
lites with Inclinations <28° 

Solutions 

GEM 3 

GEM 4 

S.E. II 

GEM 1 

GEM 2 

RMS x 10 9 

9.5 

7.6 

16.3 

22.5 

9.1 


In GEM 2 the satellite optical data (120,000obs) of GEM 1 are combined with 
gravimetry data, while in GEM 4 the satellite electronic, laser and optical data 
(400,000obs) of GEM 3 are combined with the gravimetry data. 

The agreement of GEM 3 and the French zonals with an rms value of 9. 5 
x 10' 9 for the zonal differences is quite remarkable. This is particularly so con- 
sidering that the GEM 3 zonal recovery is based upon short term zonal effects 
on weekly orbital arcs , while those of the French are based upon long term zonal 
effects. The rms of 7. 6 x 10' 9 between GEM 4 and the French zonals is ap- 
proaching the accuracy estimates as given by their standard deviations from each 
of the solutions. An rms of the zonal standard deviations from the GEM 4 solu- 
tion is 5 x 10 -9 and that for the French( 2 ), as obtained from their standard devia- 
tions, is 4. 5 x 1(T 9 for the normalized coefficients. Satellite test results using 
these solutions are presented in Section 4 for long term zonal effects. 

Geoid height zonal profiles between the French and GEM 4 solutions showed 
a maximum difference of 0. 8 meters. An estimated rms of height differences was 
less than 0. 3 of a meter. The GEM 2 and SAO S.E. II geoid height profiles, pre- 
sented in Figures 2 and 3 of the previous report, (1) show a maximum difference 
of about 1. 5 meters. 
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2. Station Height Above Geoid vs Survey and Ellipsoid Scale 


Station height above the geoid from the dynamic solution, GEM 4, is com- 
pared with the mean sea level height (MSLH) from survey and plotted in Figure 
1. Geocentric station coordinates were determined for 61 stations. The geo- 
detic station coordinates, referenced to the ellipsoid of a e = 6378155 m. and l/f= 
298.255, and geoid data are given in Table 5. As in the previous GEM 2 solution 
(Figure 6 of ref. 1), station heights from GEM 4 indicate an average equatorial 
radius for the earth of about a e = 6378145 m. In the previous results only 33 of 
the 46 stations were used, since 13 Baker-Nunn stations were excluded because 
of a problem that existed in parallactic refraction correction. This problem has 
been remedied for the present results. . 

A value of Ag 0 = 3. 5mgal, adjustment to the reference value of mean gravity, 
was obtained from use of the surface gravity data as in the previous GEM 2 solu- 
tion. The new value of Ag 0 (see the analysis on page 12 of ref. 1) corresponds to 
a scale of a e = 6378141 m. The previous result of Ag 0 = 3.3 for the GEM 2 
solution corresponded to a scale of a e = 6378142 m. The results of a e from both 
the station data and surface gravity data are reasonably consistent. Again as in 
reference 1 , the above results are based upon a fixed reference value of GM = 
398601. 3 km 3 /sec 2 . 


3. Comparison of Solutions with 5° x 5° Mean Gravity Anomalies 


A statistical procedure was presented by Kaula^ for testing satellite data 
solutions with gravimetry data and for comparing combined satellite and gravi- 
metric solutions. The following quantities are defined for the 5° x 5° mean 
gravity anomalies and used in Table 6 for comparing solutions: 


e((g t -g s ) j ) 


e(g t 2 ) 

E(G S 2 ) 

e(g t g s ) 

E(e T 2 ) 

E (V) 

E(e s 2 ) 


mean square difference between the terrestrial anomaly 
G t and the computed anomaly G s from the solution 

mean square of the terrestrial sample anomalies 

mean square of anomaly computed from the solution 

estimate of the variance of G H , the true contribution to G s 

mean square value of terrestrial anomaly error 

mean square value of neglected higher degree terms in the 
G s set (Omission error) 

mean square error in the solution G s (Commission error) 
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For a given argument Q, the preceding quantities E(Q) are computed from 


E(Q) = 


H Qi cos 
_i 

Z] COS 
i 


where the subscript i corresponds to a 5° x 5° block at latitude and cos 0j 
provides for equal area contribution over the earth. 


Formulas for the mean square errors of e T and e s and the neglected higher 
degree terms 5 g are 


E(e T 2 ) = E 


|E(G t 2 ) 


E (e s 2) = E(G S 2 ) - E(G h 2 ), 


E(5 g 2 ) = e((G t - G s ) 2 ) - E (e x 2 ) - E (e s 2 ). 


where n is the number of 1° x 1° anomalies in a 5° x 5° block. 


The gravimetry data, employed in the comparisons, consisted of 1707 blocks 
of 5° x 5° mean gravity anomalies (G T ) out of a possible 2592 world wide blocks. 
The data is listed in the appendix with a description of the reference system and 
source of data. The 5° x 5° means were based upon approximately 21,000 1° x 
1° mean free-air gravity anomalies. Each 5° x 5° mean (G x ) was formed from 
a straight average of the number, n, of 1° x 1° means that existed in the 5° x 5° 
block. The content n for the 5° x 5° blocks was quite uneven. Samples of the 
5° x 5° mean anomalies (Gx) have been selected for blocks that contain at least 
5, 10, 15, 20, and 25 points (n) of the original 1° x 1° mean gravity anomalies. 
The following table is presented for description of the 5 samples and one addi- 
tional sample containing all of the 5° x 5° means. 
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The mean square terrestrial 5° x 5° anomaly e(G t 2 ) is noted to be quite large 
in the first sample as compared to the remaining samples where e(g t 2 ) shows 
a gradual decrease. The large value (660) may be expected since the e(g t 2 ) 
for the 1° x 1° anomalies is 33 2 (~1000) mgal 2 and contain higher frequencies. 
Thus fewer points (n) used in forming a 5° x 5° mean anomaly will lack 
the effect of the ideal smoothing (i.e. with 25 1° x 1° means, n = 25) and tend to 
give large estimates for the 5 6 x 5° mean anomaly as shown in the table. In- 
spection of the 1707 blocks revealed approximately 34 5° x 5° anomalies for 
which 70 <|G t | < 219 mgal and which were formed from one or just a few points 
of the 1° x 1° original data. After these points were removed for the case of 
n> 1, e(g t 2 ) = 453 mgal 2 . The sample of n > 1 was not used in the following 
comparisons . 

The comparisons of the terrestrial 5° x 5° gravity anomaly, (G t ), with that 
computed (G s ) from different solutions are summarized in the table below. Re- 
sults for the complete set of statistics, defined by Kaula, are given in Table 6. 
The geopotential solutions used are reviewed as follows : 

GEM 1 (12 x 12) - optical satellite data solution (12 x 12 complete in 

spherical harmonics) 

GEM 3 (12 x 12) - optical plus electronic satellite data solution 

GEM 4 (16 x 16) - combined satellite (GEM 3) an4 surface gravity data 

solution 

20 x 20 - same as GEM 4 but with surface gravity complete to 

20 x 20 
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20 x 20 W2 


- same as 20 x 20 but with the surface gravity given 
twice the weight 


SAO S.E. II (16 x 16) - SAO 1969 Standard Earth II model, a combination 

of satellite optical data and surface gravity data 
consisting of 935 (equal area 300 nm square) blocks 
of mean terrestrial gravity anomalies. 

The above solutions all contain selected higher degree zonal and satellite 
resonant coefficients extending to degree 22. 

Comparison of 5° x 5° Mean Gravity Anomalies Obtained from Potential 
Coefficients (G s ) and Terrestrial Data (G T ) (Summary of Table 6) 


Mean Square of Differences (G x - G s ), E^(g t - G s ) 2 ^ mgals 2 


Solutions 

n>5 

1284 (Blocks) 

n>10 

881 

n>15 

624 

n>20 

434 

n = 25 
233 

GEM 1 (12 x 12) 

261 

242 

• 249 

216 

213 

SAO S.E. II (16 x 16) 

261 

241 

242 

200 

224 

GEM 3 (12 x 12) 

252 

231 


208 

204 

GEM 4 (16 x 16) 

208 

186 


164 

160 

20 x 20 

179 

153 


130 

125 

20 x 20 W2 

172 

147 


123 

119 

e(g t 2 ) 

435 H 



367 

374 

E(e T 2 ) 

55 



23 

23 


Variance of the Error e s of Commission in the Solution, E (e s 2 ) mgal 2 



n>5 

n>10 

n>15 

n>20 

n = 25 

GEM 1 (12 x 12) 

24 

23 

26 

30 

22 

SAO S.E. II (16 x 16) 

36 

39 

43 

38 

41 

GEM 3 (12 x 12) v „ 

21 

18 

21 

26 

18 

GEM 4 (16 x 16) ’I-*’ 

■ V 

■ 14 


6 

12 

: i 

4 


7 

























Variance of Error of Omission 5 g for the Solution, E(5 g 2 ) mgal 2 



n>5 

n>10 

n>15 

n >20 

n = 25 

GEM 3 (12 x 12) 

176 

177 

190 

159 

163 

GEM 4 (16 x 16) 

137 

143 

154 

130 

133 

20 x 20 W2 

104 

106 

120 

101 

94 


The statistical results in the table are in general quite good. In each of 
the solutions, the sample variances e((G t - G s ) 2 ) generally show better results 
for the samples where the 5° x 5° mean gravity anomalies are based upon a more 
complete set of points (n) of the 1° x 1° means . This may be expected for the 
solutions where the surface gravity data is used, since the 5° x 5° mean anom- 
alies are weighted proportionately to the number of points n (n + 1 is actually 
used). The satellite solutions, however, also show the same trend of agreement. 
The variances show the best agreement for the 20 x 20 field W2, where the sur- 
face gravity data is given greater weight by a factor of 2. The increased weight- 
ing contributes a reduction of 6 to 7 mgal 2 for the associated samples w.r.t. 

20 x 20 field which has the normal weight. However the increase in the dimen- 
sion to degree and order 20 provides a larger reduction of 29 to 36 mgal 2 w.r.t. 
the GEM 4 field (16 x 16). And in turn it has a reduction on the average of about 
50mgal 2 over GEM 3 (12 x 12) satellite solution. The GEM 3 satellite solution 
(electronic plus optical data) is seen to agree better with the gravity data than 
GEM 1 (optical data only) satellite solution in every sample with a reduction in 
the variances of from 8 to 11 mgal 2 . 

The SAO S.E. II solution, based upon satellite optical and surface gravity 
data of 300 nm equal area squares, have variances which are a little larger 
than the GEM 3 satellite solution. This result is not expected, since surface 
gravity data was used in the S. E. II and not used in the GEM 3 satellite solution. 

An estimate for the accuracy of the solution is given in the table by the vari- 
ances E(e s 2 ). These error estimates presented for the GEM 1, GEM 3, and 
GEM 4 solutions decrease respectively as may be expected from the above agree- 
ment with the surface gravity data as indicated by e((Gt - Gs) 2 ). The error es- 
timates for the S.E. II solution are larger than the other solutions. The error 
estimates for each solution are relatively consistent on each of the 5 samples. 

Since 5° x 5° mean gravity anomalies covering the earth contain information 
that may be expected to correspond td a geopotential solution with spherical har- 
monics complete to degree and order 36, the neglected portion 5 g of the solu- 
tions should be expected to contain much remaining information. The information 
that remains to be recovered in the data as given by the variances, E (5 g 2 ), is 
relatively consistent on each of the samples for the GEM 3, GEM 4, and the 
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20 x 20 W2 solutions which are presented in the table. The sample variances 
for each of the solutions , where 

Total _ Data Solution Omission 
Mis closure Error Error Error 

e((G t - G s )2) - E(e T 2) + E(e s 2) + E (5 g 2 ) 

show that the larger part of the misclosure lies in the omitted part of the ex- 
pected solution of the terrestrial data (36 x 36), namely E(5 g 2 ). The 20 x 20 W2 
field shows approximately lOOmgal 2 that remains to be recovered from the data, 
while the GEM 3 satellite solution shows this value as large as 175mgal 2 and 
GEM 4 solution shows 137 mgal 2 . 

The degree variances of the gravity anomalies are presented in Table 7 for 
the different solutions (for GEM 1 see Table 6 of ref. 1). The sum of the degree 
variances of the gravity anomalies should be equal to e(g s 2 ) over the entire 
Earth. Values of e(g s 2 ) are given below for the case of n>5 (1284 blocks) and 
for 2592 5° x 5° blocks covering the entire earth. 


Solution 

e(g s 2 ) 

n>5 1284 Blocks 

e(g s 2 ) 

2592 Blocks 

Sum of Degree 
Variances 
(n = 2 to 22) 

GEM 1 (12 x 12) 

222 

186 

186 

GEM 3 (12 x 12) 

225 

191 

188 

GEM 4 (16 x 16) 

257 

240 

234 

SAO S.E. II (16 x 16) 

248 

219 

1 

216 


The results of e(g s 2 ) for the entire earth show good agreement with the sum of 
the degree variances. The values of e(g s 2 ) for the sample anomalies correspond- 
ing to the case of n>5 all are larger than the other two cases. Since this is true 
for the satellite solutions as well, it indicates that the areas not sampled, prin- 
cipally in the South Pacific , have smaller gravity anomalies at least as repre- 
sented by the spherical harmonic solutions. 

4. Satellite Test Results 

Two geopotential solutions have been tested on all 23 of the satellites (Table 
3A) with optical data. Results of the rms of observation residuals are listed in 
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Table 8 for a weekly arc on each of the satellites. Similar rms values, given 
in Table 3A, are based upon the starting solution. Solutions for which the results 
were obtained are the GEM 1 and the S.E. II models. The average rms for all 
the satellites for each model is as follows: 

SAP S.E. II GEM 1 

Average rms 3.42" 2.75" 

These results are somewhat expected since the SAO S.E. II solution also had to 
satisfy surface gravity data. On the other hand, the GEM 1 satellite solution 
satisfied the GSFC surface gravity data almost as well as the SAO S.E. II solu- 
tion (see Table 6). 

A result is given below for the GEM 3 and SAO S.E. II solutions employing 
ISAGEX French laser data of 1700 range observations, taken on the Haute 
Provence site for a weekly arc of the GEOS-I satellite (710219/26). This data 
was not used in either solution, and the rms of the laser range residuals are: 

S.E. II GEM 3 


rms 8 meters 4 meters 1700 range obs 

Zonal solutions were tested by Carl A. Wagner^) on 21 satellites, including 
the low inclination satellites of SAS and PEOLE , for their long term zonal ef- 
fects on mean elements. Wagner uses as a test criteria the weighted rms of the 
mean element residuals of each solution. These residuals are to be used as a 
relative measure of testing and are listed here with other solutions (not included 
in ref. 4) for comparison as follows: 


Solution 

rms 

SAO S.E. II 

5.49 

GEM 2 

4.80 

GEM 1 

3.62 

French 71 

3.28 

GEM 3 

2.92 

GEM 4 

2.89 

Wagner 

1.50 


Only about half of the 21 satellites used in Wagner's solution are contained in 
the other models. 
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The GEM 3 satellite solution and the GEM 4 satellite/gravimetric solution 
compare well in this test. Considering the GEM 3 solution was based upon 
satellite weekly orbital arcs and the French 71 solution was based upon satellite 
long term zonal effects, it verifies that good zonal recovery may be achieved 
from short term zonal effects. The SAO S.E. II, GEM 1 and GEM 2 solutions 
do not compare as well in these tests because they do not contain the effects of 
low inclination satellites. Wagner’s result is expected to have the lowest rms 
since his solution is based entirely upon the test data. 
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III. SUMMARY AND CONCLUSIONS 

The new satellite solution GEM 3 including the addition of two low inclination 
satellites and employing some 400,000 observations of electronic, laser and 
optical data provides for a refinement over the GEM 1 satellite solution. GEM 
1 was based upon 25 satellites consisting primarily of optical data. The test 
results of GEM 3 with gravimetric data and long term satellite zonal effects show 
improvement over the GEM 1 solution and also that of the SAO 1969 S.E. II 
model. The latter two models contained satellites all with inclinations greater 
than 28°. A test with GEOS-I ISAGEX French laser data, independent of the 
solutions, also showed better results for the GEM 3 than the S.E. II model. 

The GEM 4 solution combined the GEM 3 satellite data with the gravimetric 
data and provided a geopotential model complete to degree and order 16 with 
zonals and selected satellite resonant coefficients extending to degree 22. The 
solution included 61 center of mass tracking station locations and an adjustment 
(Ag 0 ) to the reference value of mean gravity. Analysis of the heights of these 
stations above the GEM 4 geoid with the mean sea level heights from survey in- 
dicated a mean earth ellipsoid radius of about a e = 6378145 meters, while analy- 
sis of Ag 0 (3 . 5 mgal) indicated a e = 6378141 meters. The two results are fairly 
consistent and are based upon the reference value of GM = 398601. 3km 3 /sec 2 . 

The analysis of the gravity data in terms of 5° x 5° mean gravity anomalies 
indicated that additional information remains to be recovered from the data. A 
GSFC combination satellite/gravimetric solution complete to degree and order 
20 showed approximately 100 mgal 2 remained to be recovered for the 5° x 5° 
mean terrestrial anomaly. 
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Table 1 


Comparison of Goddard Earth Models (GEM) 
(Geopotential Solutions) 


Geopotential 

Solution 

Completeness of 
Spher. Harm. 
Field* 

Satellite Data 

Gravimetry Data 

GEM 1 

12 x 12 
Satellite 

120, 000 Optical Obs. 
on 25 Satellites 


GEM 2 

16 x 16 
Combined 

GEM 1 Data 

1707 5° x 5° Mean 
Gravity Anomalies 

GEM 3 

12 x 12 
Satellite 

400, 000 Optical and 
Electronic Obs. on 
27 Satellites 


GEM 4 

16 x 16 
Combined 

GEM 3 Data 

1707 5° x 5° Mean 
Gravity Anomalies 

*A11 solutions includ 

ed higher degree zonal and satellite resonant coefficients extending to degree 22 . 


(GEM 1 and 2 excluded the zonal coefficient of degree 22.) 

• 61 Center of mass tracking station locations included in the GEM 
3 and 4 models . 

46 stations included in the GEM 1 and 2 models 

• Two low inclination satellites SAS 3°, PEOLE 15° included in 
GEM 3 and 4. 

25 satellites in previous GEM 1 and 2 solutions all have inclination 
>28°. 

• 5° x 5° mean gravity anomalies based upon: 

19,000 1° x 1° Anomalies from ACIC 

2,000 1° x 1° Anomalies from other sources 
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GSFC Geopotential Solutions (Normalized Coefficients x 10 6 ) 1 


O a* n n 

Q O (MO 

♦ cm mm 

u cy n - 


« S S •C 

9 o 9 e 

M « CD M 

or* o cm 


w in - 4 a 

® - n cm v 9 

o> - o m <o o 

cm in on: -JO 


0> CO ® B 

O' <# ~ •* 

no ^ in 

O C JO 


n in o « 

ns -on 

<m - - n 

o u o o 


••■> 4 O X - J 

n <o s * in *■ 

(M >0 9 V -0 CM 

n - O IM wN 


win in m w in in in ww ww win tnw win 


m in -o <j -o 'V u -o <o« 


uw y ® *> ® u ® u m u ® u in y ® y ® u ® y ® y ® y ® uw y ® 


n <t> o><o - r> cOi 

_ O' - OO -O' *M W 

” 9 in « cj - — ® 

5 — — o — o— — o 


w n n n 

— CM — W 

OO - CM CM 

it- - n « w 


• O' — n -o <d -* -* — o in — 

i a no* cm n oo o o>s or* 

• a in-o c n o'-* - s — — 

1 w (MO OO C — O O CO- 


WS O' O' 

9 ® ® |\J 

cu <o n -a 

O O - O 


OO 0 J OO OO 


s - in « 

J HI JO 

- n „ w 

no - <* 


O' ■» j ® nj O' -at 

n ® on - n nr 

O' O 9 «■ <\l - O' . 

CM — CM O O O O ■ 


ss n a ns 09 

n n cm o <t n 

oo <* — <ns no 

oo oo oo ou 


u ® u ® y ® u ® 


> ® y ® y ® y ® o® 


cm s o ® in s 

9 r> ^ ® s ® 

w <o in — ® o 

n ® an on 


CM O ® -0 m ® 

® W ns so 

9 0 cm - <* <a 

O O C’ O © o 


on o s 

n -* -« n 

o - o <•> 


ns in a 

c O' in d 

CM - s -* 

S 9 O' CM 


i o a o w 9 ® - 

■ s cm w w cm n n 

) S CM CM a CM CM ® 

C CM w O CM CM CO 


> o ® n. n - ns - o >s ® - 

> — o- an — ® OO a cm w a 

. w — CM U <M fv-. S® — W CM 

> — J ■> S-sf 9 CM 9 CM OO (MM 


nn nn nn nn 


win a a nn <r <a in w a a 


i ® u ® o® u ® o® u® u® o® 


— ® T ® O CM 

o o in cm w w 

s — n ® — 

— c t o c 


9 0 9 — -o W-* n<* S® ®® »"-» 

-o® n— -tf® o ® «n n c?s no wn 

— <* no ss o>o w n cmw w>> s <t 

O CM W 9 O •-* O •-* CM — JO - J O - 

(Mo JO O O a *.J O J J '5 O O JO 


o— >m® nn n— no <»<m s® *cm 

o»s o>® aw ®c- <»— — ® no* ®cm 

«n w<r — s c n c-» no cm® — n 

jo o o jo o o — -j o j <» n o « 


in® <oo s s to® 


nn <»<» u*W oo cm cm 


u® o® o® u® o® o® u® u® u® u® o® y® u® o® u® u® u® u® u® 


OC Oo C O vJO OO OO CO 


>o CO OC OO OO OO OO OO OO OO 


> 'J 0 9 JO O O 


in 9 9 ® ® c <J 


U® U® u® u® U® O ® y® y® y® o® y® y® y® y® y® u ® y® y® u® 


16 



GSFC Geopotential Solutions (Normalized Coefficients x 10 6 ) 2 
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Satellite Optical Data 

300 Weekly Satellite Arcs Opt. Data (Primarily SAO Baker-Nunn) 
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Satellite Electronic, Laser, and Additional Optical Data (48 Weekly Orbital Arcs) 
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^Observations (passes) 


Table 4 


Comparison of Zonal Coefficients (Normalized Coefficients x 10 6 ) 


REE 

FRENCH 

1971 

GEM 3 

GEM 4 

SAO 

1969 

GEM 1 

GEM 2 

2 

- 484.170 

- 484.171 

- 484.169 

- 484.166 

- 484.177 

- 484.167 

3 

.961 

.958 

.957 

.959 

.962 

.955 

4 

.540 

.547 

.541 

.531 

.557 

.537 

5 

.068 

.068 

.069 

.069 

.062 

.073 

6 

-.155 

-.162 

-.153 

-.139 

-.178 

-.145 

7 

.094 

.092 

.091 

.094 

.105 

.087 

8 

.051 

.062 

.051 

.029 

.080 

.040 

9 

.027 

.030 

.031 

.023 

.008 

.033 

10 

.051 

.040 

.050 

.077 

.021 

.065 

11 

-.049 

-.056 

-.056 

-.042 

-.020 

-.055 

12 

.038 

.046 

.039 

.008 

.059 

.021 

13 

.039 

.049 

.048 

.024 

.002 

.043 

14 

-.016 

-.030 

-.027 

.014 

-.037 

-.009 

15 

.015 

-.007 

-.005 

.031 

.047 

.004 

16 

-.008 

-.012 

-.009 

-.033 

-.013 

-.026 

17 

.005 

.020 

.017 

.014 

-.035 

.007 

18 

.023 

.016 

.011 

.038 

.018 

.023 

19 

.018 

.008 

.009 

.035 

.045 

.015 

20 

.014 

.003 

.009 

.001 

-.002 

-.001 

21 

-.016 

-.008 

-.008 

-.022 

-.031 

-.012 

22 


-.001 

-.004 





Differences with French Zonals x 10 9 


GREE 

GEM 3 

GEM 4 

SAO 

GEM 1 

GEM 2 

2 

- 1 

1 

4 

- 7 

3 

3 

- 3 

- 4 

- 2 

1 

- 6 

4 

7 

1 

- 9 

17 

- 3 

5 

0 

1 

1 

- 6 

5 

6 

- 7 

2 

16 

-23 

10 

7 

- 2 

- 3 

0 

11 

- 7 

8 

11 

0 

22 

29 

-11 

9 

3 

4 

- 4 

-19 

6 

10 

-11 

1 

26 

-30 

14 

11 

7 

- 7 

7 

29 

- 6 

12 

8 

1 

-30 

21 

-17 

13 

10 

9 

-15 

-37 

4 

14 

-14 

-11 

30 

-21 

7 

15 

-22 

-20 

16 

32 

-11 

16 

- 4 

- 1 

-25 

- 5 

18 

17 

15 

12 

9 

40 

2 

18 

- 7 

-12 

15 

- 5 

0 

19 

-10 

- 9 

17 

27 

3 

20 

-11 

- 5 

•13 

-16 

15 

21 

8 

8 

- 6 

-15 

4 

RMS 

9.5 

7.6 

16.3 

22.5 

9.1 
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Table 5 


GEM 4 Station Coordinate Solution 


GEOIO 


TYPE 

STATION 

L ATITUOF 

LONGI TODE 

HEIGHT 

MSL 

HEIGHT 

RES 


NAME NUMBER 

DEG MIN SEC 

DEG MIN SEC 

M 

M 

M 

M 


M 

1HP0IN 

1021 

38 

25 

49.920 

282 

54 

48.543 

-42.3 

5.8 

-42.3 

-5*8 

M 

IFTMYR 

1022 

26 

32 

53.408 

278 

8 

4.219 

-32.3 

4.8 

-29.0 

-7.3 

M 

IPGMER 

1024 

-31 

23 

25.122 

136 

52 

15.442 

129.5 

132.8 

-3.0 

-0.3 

M 

1SATAG 

1020 

-33 

0 

58.371 

289 

19 

53.674 

707.9 

693.4 

23.7 

-9.2 

M 

1M0JAV 

1030 

35 

19 

47.898 

243 

5 

59.184 

889.5 

929.1 

-28.9 

-10.7 

M 

1JOBUR 

1031 

-25 

53 

0.742 

27 

42 

26.375 

1537.3 

1522.0 

22.3 

-7.0 

M 

INEWFL 

103? 

47 

44 

29.861 

307 

16 

46.193 

63.2 

69.0 

12.0 

-17.8 

M 

1GFORK 

1034 

48 

1 

21.333 

262 

59 

19.474 

215.2 

252.6 

-30.1 

-7.3 

M 

1WNKFL 

1035 

51 

26 

46.074 

359 

18 

8. 307 

104.4 

67.4 

45.3 

-8.3 

M 

IROSMN 

1037 

35 

12 

7.395 

277 

7 

41.380 

867.3 

909.3 

-37.9 

-4.1 

M 

IORORL 

1038 

-35 

37 

32.109 

148 

57 

14.881 

944. 

931.6 

15.0 

-3.2 

M 

IROSMA 

1042 

35 

12 

7.398 

277 

7 

41.169 

867.5 

909.4 

-37.9 

-4.0 

M 

1 TAN AN 

1043 

-19 

0 

31.629 

47 

17 

59.444 

1362.6 

1378.0 

-S.l 

-9.3 

M 

1UNDAK 

7034 

48 

1 

21.403 

26? 

59 

19.464 

216.7 

252.6 

-30.1 

-5.8 

M 

1 EDI NB 

7036 

26 

22 

46.796 

261 

40 

7.504 

21.3 

59.6 

-18.9 

-19.4 

M 

I COL BA 

7037 

38 

53 

36.235 

267 

47 

40.934 

227.6 

272.7 

-32.3 

-12.8 

M 

1BERMD 

7039 

32 

21 

49.897 

295 

20 

35.029 

-17.1 

31.2 

-43.4 

-4.9 

M 

1 PUR 10 

7040 

18 

15 

28.918 

2 94 

0 

23.629 

-7.4 

49.7 

-47.5 

-9.6 

M 

1GSFCP 

7043 

39 

1 

15.736 

263 

10 

20.326 

3.0 

53.5 

-41.9 

-8.6 

M 

1DENVR 

7045 

39 

38 

48.050 

255 

23 

38.668 

1758.2 

1790.0 

-21.5 

-10.3 

M 

IJUM40 

7072 

27 

1 

14.408 

279 

53 

12.876 

-33.6 

14.2 

-34.2 

-13.6 

M 

1SUDBR 

7075 

46 

27 

21.244 

279 

3 

10.410 

236.0 

281.9 

-39.0 

-6.9 

M 

1J4MAC 

7076 

18 

4 

34.762 

283 

11 

27.073 

419.6 

445.9 

-25.2 

-l.i 

B 

10RGAN 

9001 

32 

25 

24.986 

253 

26 

49.004 

1617.3 

1651.0 

-23.9 

-9.8 

B 

10LFAN 

9002 

-25 

57 

36.038 

28 

14 

52.527 

1561.0 

1544.0 

21.6 

-4.6 

B 

100MER 

9003 

-31 

6 

2. 109 

136 

47 

3.350 

153.7 

162.5 

-2.7 

-6.1 

B 

1 SPA I N 

9004 

36 

27 

46.734 

353 

47 

36.956 

57.3 

25.9 

50.0 

-18.6 

B 

l TOKYO 

9005 

35 

40 

22.975 

139 

32 

16.555 

83.9 

59.8 

38.2 

-14.1 

B 

INATOL 

9006 

29 

21 

34.687 

79 

27 

27.520 

1869.4 

1927.0 

-50.2 

-7.4 

B 

IOUIPA 

9007 

-16 

27 

56.824 

288 

30 

24.601 

2484.7 

2452.0 

29.6 

3.1 

B 

1SHRAZ 

9008 

29 

38 

13.743 

52 

31 

11.369 

1578.7 

1596.0 

-9.0 

-8.3 

B 

l CUR AC 

9009 

12 

5 

25.056 

291 

9 

44.547 

-24.4 

8.7 

-29.3 

-3.8 

B 

UUP TR 

9010 

27 

1 

14.020 

279 

53 

13.375 

-24.9 

15.1 

-34.2 

-5.8 

B 

1 V ILDO 

9011 

-31 

56 

34.718 

294 

53 

36.624 

626.6 

598.4 

24.6 

3.6 

B 

IMAUIO 

9012 

20 

42 

26.051 

203 

44 

33.984 

3040.4 

3034.0 

4.5 

i. 9 

B 

HOPKIN 

9021 

31 

41 

3.174 

249 

7 

18.592 

2337.2 

2382.0 

-26.7 

-18.1 

B 

AUSBAK 

9023 

-31 

23 

25.874 

136 

52 

43.675 

136.6 

141.2 

-3.0 

-1.6 

B 

DEZEIT 

9028 

8 

44 

51.106 

38 

57 

33.366 

1900.8 

1924.0 

-6.7 

-16.5 

B 

COMRIV 

9031 

-45 

53 

12.418 

292 

23 

9.422 

195.3 

186.5 

11.2 

-2.4 

* B 

GREECE 

9091 

38 

4 

44.736 

23 

55 

58.650 

483.2 

467.0 

32.1 

-15.9 

♦ B 

EDWAFB 

9425 

34 

57 

50.520 

242 

5 

7.955 

747.3 

784.2 

-30.9 

-6.0 

* B 

JOHNST 

9427 

16 

44 

38.733 

190 

29 

9.342 

17.3 

5.0 

10.9 

1.4 

G 

MADGAR 

1123 

-19 

1 

14.228 

47 

18 

11.362 

1386.9 

1399.0 

-6.1 

-6.0 

G 

ROSRAN 

1126 

35 

11 

45.641 

277 

7 

26.252 

827.0 

873.9 

-37.9 

-9.0 

G 

ULASKR 

1128 

64 

58 

19.069 

21 2 

29 

12.706 

341.4 

346.6 

10.5 

-15.7 

G 

CARVON 

1152 

-24 

54 

10.766 

113 

42 

59.736 

6.1 

37.9 

-20.6 

-11.2 

L 

GOOLAS 

7050 

39 

1 

14.539 

283 

10 

18.706 

13.4 

54.8 

-41.9 

0.5 

L 

WALLAS 

7052 

37 

51 

36.235 

284 

29 

23.889 

-42.9 

8.6 

-43.1 

-8.4 

L 

CRMLAS 

7054 

-24 

54 

15.794 

113 

42 

58.262 

-6.6 

31.4 

-20.6, 

-17.4 

♦ D 

ANCHOR 

2014 

61 

17 

0.152 

210 

10 

28.960 

64.0 

68.0 

15.0 

-19.0 

D 

TAFUNA 

2017 

-14 

19 

50.005 

189 

17 

3.159 

27.5 

6.1 

25.5 

-5.1 

* D 

WAHIWA 

2100 

21 

31 

15.577 

202 

0 

10.406 

396.2 

388.0 

7.6 

0.6 

♦ D 

L ACRES 

2103 

32 

16 

44.546 

253 

14 

45.439 

1151.6 

1203.0 

-24.0 

-27.4 

* 0 

LASHAM 

2106 

51 

11 

9.350 

358 

58 

25.632 

218.5 

190.3 

45.7 

-17.5 

♦ 0 

APLMND 

2111 

39 

9 

48.630 

283 

6 

11.907 

95.7 

145.0 

-41.8 

-7.5 

* D 

THOLEG 

2018 

76 

32 

20.076 

291 

13 

52.887 

53.7 

43.0 

13.8 

-6.1 

♦ D 

PRETOR 

2115 

-25 

56 

48.233 

28 

20 

52.072 

1581.0 

1580.0 

21.4 

-20.4 

♦ D 

ASAMOA 

2117 

-14 

19 

50.098 

189 

17 

2.960 

38.5 

6.0 

26.5 

6.0 

D 

MESHED 

2817 

36 

14 

26.605 

59 

37 

44.419 

953.1 

994.6 

-19.9 

-21.6 

0 

FRTLMY 

2822 

12 

7 

54.040 

15 

2 

6.939 

299.1 

298.4 

13.0 

-12.3 

D 

NBRZIL 

2837 

-5 

54 

57.855 

324 

49 

55.940 

6.2 

41.0 

-11.2 

-23.6 


'MEAN SEA LEVEL HEIGHTS. MSL, FOR THESE SITES MAY NOT BE RELIABLE. 
M-MOTS. B-BAKER-NUNN, G-GRARR, L-LASER, D- NWL DOPPLER. 
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Comparison of Satellite and Combination Solutions with 
Surface Gravity Measurements (mgal 2 ) 
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Table 6 (continued) 
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Table 7 


a n 2 Degree Variances of Gravity Anomalies 


Degree n 

GEM 3 

GEM 4 

SAO 69 

20 x 20 
Field 

0 

— 

3.5 

2. i 

3.0 

2 

7.3 

7.3 

7.3 

7.3 

3 

33.707 

33.699 

32.842 

33.717 

4 

21.317 

21.343 

21.805 

21.312 

5 

21.895 

21.968 

17.785 

21.965 

6 

19.029 

19.324 

15.652 

19.263 

7 


19.287 

15.491 

19.211 

8 


10.117 

6.639 

10.081 

9 


11.347 


11.406 

10 


10.438 

12.860 

10.473 

11 

7.902 

8.375 

12.234 

8.059 

12 

5.201 

5.360 

5.099 

5.595 

13 

4.222 

17.586 

11. 121 

19.436 

14 i 

1.739 

14.818 


14.986 

15 

1.097 

17.230 


17.380 

16 

0.600 

6.911 

13.844 

7.078 

17 

1.081 

0. 995 

2.253 

19.307 

18 

3.432 

3.152 

1.667 

18.166 

19 

0.373 

0.409 

3.545 

25.274 

20 

2.186 

1.870 

1.693 

9.646 

21 

0.957 

0.968 

0.185 

0.763 

22 

2.620 

2.496 

0.0 

2.076 










































Table 8 


Weighted RMS of Optical Observation Residuals on 23 Satellites 
(weight = 1/2 of a second of arc) 


Satellites 

SAOS.E. II 

GEM 1 

TELSTAR 

2.46 

0.95 

GEOS-A 

1.08 

0.92 

SECOR 

1.36 

1.27 

OV1-2 

2.10 

2.06 

ECHO 

1.31 

1.16 

Dl-D 

2.50 

1.66 

BE-C 

1.16 

0.92 

Dl-C 

1.88 

1.08 

ANNA -IB 

1.33 

1.16 

GEOS-B 

1.22 

0.90 

OSCAR 

1.34 

1.17 

5BN-2 

2.61 

2.56 

COURIER 

1.28 

1.44 ! 

GRS 

2.64 

1.89 

TRANSIT 

1.33 

1.23 

BE-B 

1.44 

1.29 

OGO-2 

2.93 

1.49 

INJUN 

2.00 

2.28 

AGENA 

2.45 

1.88 

MIDAS 

0.87 

0.86 j 

VANG-2R 

0.87 

0.80 

VANG-2S 

1.53 

1.42 

VANG-3S 

1.64 

1.26 

TOTAL 

39.33 

31.65 

Average mas 
(in seconds of arc) 

3.42" 

2.75" 
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APPENDIX 


DESCRIPTION AND TABULATION OF 5° x 5° MEAN GRAVITY ANOMALIES 

The source of most of the gravimetry data was the U.S. Aeronautical Chart 
and Information Center (5) which provided 19,000 one-degree by one-degree 
mean free-air gravity anomalies. A further set of 2000 mean gravity anomalies 
were obtained from a number of other sources. These data were used to form 
1707 five-degree by five-degree mean gravity anomalies by a straight averaging 
of the one-degree by one-degree mean anomalies , and provided a total coverage 
of about 70% of the earth's surface. (See Figure 1 in reference 1 for a map of 
surface gravity data coverage.) 

Table A-l lists the 1707 5°-by-5° mean anomalies represented at the mid- 
point of the 5°-by-5° block in latitude and longitude. The mean anomalies listed 
are referred to the ellipsoid and normal gravity whose parameters are: 

a e = 6378155 meters 

f = 1/298.255 

GM = 3.986013 x 10 14 m 3 /sec 2 

a) = 0. 7292115146 x 10 -4 radians/sec 

The original gravity anomalies Ag It referred to the International Gravity 
Formula ( 7 j), were converted to Ag R inthe above reference system as follows: 

go 1 = Ag! + 7 1 , 

go = go 1 + Potsdam correction, 
and Ag R = g 0 “ 7r » 

or Ag R = Agj + 7 1 “ 7 R + Potsdam correction, 

where 

7 R = 7e (1 + P sin 2 0 + 0* sin 2 20). 


A-l 



The above quantities are defined as follows : 


go 1 denotes measured gravity reduced to the geoid (mean sea-level) in the 
Potsdam system, 

Ag\ is the gravity anomaly in the International system, 

7i is normal gravity on the International Ellipsoid, 

g 0 is measured gravity reduced to the geoid in the absolute system, 

Potsdam correction is the constant correction needed to convert Potsdam 
system values to an absolute system. 


7r is normal gravity on the reference ellipsoid, and 
Ag R is the mean anomaly in the reference system. 

The equatorial gravity (7 e ) and the constants 0 and 0* are obtained from the 
reference parameter of GM, a e , f, andcoby the following relations:^) 


GM = a 2 (l-f)7 e ^l+-|m+|-fm+-| m 2 j, 
0 = Pi + 02 

?* =-j& 


Hi 


= -f + | m + 



7 4 


2 


02 



m 


co 2 a 3 (1 - f) 
GM 


From the above, 

7r = 978.0291 (1.0 + 0.0053025 sin 2 # - 0.00000585 sin 2 2#) gal 


A-2 



and 


7i = 978.049 (1.0 + 0.0052884 sin 2 0 - 0. 0000059 sin 2 20) gal. 

With an adopted value for the Potsdam correction of -13.7 mgals, the con- 
verted anomalies (in mgals) become, 

Ag R = Agj + 6.2 - 13.7 sin 2 0. 

In Table A-l that follows, the columns are labeled PHI, LAMDA, DEL G, 
and N, where 

o 

PHI is the latitude of the midpoint of the 5° x 5° block, 

LAMDA is the longitude of the midpiont of the block , 

DEL G is the mean anomaly (Ag R ) in the reference system given above, and 
N is the number of 1° x 1° mean anomalies within a 5° x 5° block. 
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Table A-l 

x 5° Mean Gravity Anomalies 
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